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1. Introduction 

In this article we write the equations of barotropic compressible fluid mechanics 
as a geodesic equation on an infinite-dimensional manifold. The equations are given 

by 

(1) ut + VuU = gradp, 

P 

(2) pt + div {pu) = 0, 

where the fluid fills up a compact manifold M, m is a time-dependent velocity field 
on M, and p is the density, a positive function on M. The barotropic assumption is 
that the pressure p is some given function of the density, although our methods also 
extend to certain more general isentropic flows. Our infinite-dimensional manifold 
is the product 'D{M) x C°°(Af, R). This is a group using the semidirect product 
(which is sometimes incorporated in other treatments) , but the Ricmannian metric 
we use is neither left- nor right-invariant. Hence our geodesic equation is not an 
Euler- Arnold equation. Wc compute the sectional curvature and show that at least 
when AI ~ S^, the curvature is always nonnegative. We also establish some results 
on the Lagrangian linear stability of solutions of this system, for certain nonstcady 
solutions in one dimension and steady solutions in two dimensions. 

It has been known for many years that the system ([T|)-([2|) can be derived via 
Hamilton's principle on an infinite-dimensional manifold (see e.g., Ebin |Eb| . Smo- 
lentsev jSmj , and Holm et al. |HMR| , but all such approaches use a nonzero potential 
energy, so that the equations are essentially Newton's equation ^ = — grad ^' for 
some nonzero ^' rather than a pure geodesic motion (with = 0). We review this 
approach in Section In contrast we obtain a genuine geodesic equation, but the 
price we pay is that our geodesic equation only gives the barotropic equations on 
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a certain iionliolonomic distribution (i.e., some geodesies correspond to barotropic 
fluids, while others have no physical meaning). Our method is somewhat closer to 
the metric introduced in Eisenhart [Ei] on the product I?(M) x R^, where extra 
degrees of freedom are introduced and weighted in the Riemannian metric by terms 
involving the potential energy and an arbitrary parameter. The difference is that 
our metric is "diagonal" on the product, and the extra variables we introduce are 
directly related to the density. 

One motivation for doing this is to obtain some intuition for the behavior of com- 
pressible fluids under a perturbation. In much the same way that a one-dimensional 
particle trajectory in a convex potential energy behaves approximately like a geo- 
desic in a surface of positive curvature, we hope to understand stability of a com- 
pressible fluid using curvature computations which at least intuitively suggest the 
behavior of perturbed solutions. Another motivation is to understand the warped 
product geometry of 'D{M) x C°°(M, R), under a noninvariant metric. Invariant 
metrics have some nice algebraic properties, but are often not physically relevant 
(for example, the right-invariant metric on T>{M) does not come from any kind 
of physics, and its geodesic equation has no known physical relevance; while the 
non-invariant metric we consider is naturally related to the kinetic energy, and 
its geodesic equation describes the motion of a force- free family of particles). 

In Section [31 we write down the Riemannian metric and its geodesies along with 
its Riemannian curvature. The advantage of the geodesic approach, as pioneered by 
Arnold for incompressible fluids, is that one can in principle use the curvature 
to discuss Lagrangian stability of the fluid: intuitively, if the sectional curvature is 
positive, then the particle paths should be stable under small perturbations of the 
initial velocity. (More precisely one uses the Ranch theorem which gives results up 
to the first conjugate point.) This cannot be applied rigorously in our case, since 
unlike in the case of incompressible fiuids |EM| , the Riemannian exponential map 
cannot be smooth in any Hilbert topology: simple examples show that geodesies 
in the present context are not even locally minimizing. Hence our approach is 
necessarily only formal, and so we may as well require all objects to be C°° rather 
than working in Sobolev spaces. 

Nonetheless we can still analyze the Jacobi equation for linear perturbations of 
geodesies directly, which we do in Section H] see |BG| for a similar perspective. 
We establish stability or weak instability in special cases (for arbitrary solutions in 
one space dimension and a rigid rotational fiow in two space dimensions). This is 
not a significant drawback, since even for incompressible fiuids a direct approach 
to stability is frequently more informative than an analysis of the curvature (see 
jP2| ). In the one-dimensional case we show that Lagrangian perturbations grow 
at most linearly in time up to the shock (when solutions cease to be smooth and 
our methods no longer apply), while in the two-dimensional steady case, many 
Lagrangian perturbations are bounded for all time. This portion of the paper 
previously appeared in the author's thesis jPlj . 

2. Background 

We describe a compressible fiuid as follows: consider a manifold N with a volume 
form V (describing the mass distribution) and a Riemannian manifold M with 
metric (•,•) and Riemannian volume form ^ (describing physical space). Fluid 
configurations are described by trajectories 'q{t) G C°°{N,M), where the density is 
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defined by 

(3) porjrfii~v, or ]ac{r/)p o rj. 

This space is formalljQ a manifold, where the tangent spaces are 

Tr,C°°{N, M)^{U: N ^ TM \ U{p) e T,-,(p)M\lp e N} 

and the coordinate charts are given by the exponential maps 

Exp,,(C/) =p^ exp,,(p)(f/(p)), 

where exp is the Riemannian exponential map on M (which takes the velocity 
vector U{p) to the point on the geodesic through 77 (p) in direction U{p) at time 
one). In the Riemannian metric 



(4) ((f/,nr, - / (C/(p),F(p)),(p)di/, 

J N 

Exp is precisely the Riemannian exponential map. On the space C°°{N, M), the 
map Exp is globally defined, although an obvious problem is that physically maps 
from N to M should be one-to-one: the physically natural space is the space of 
smooth embeddings S{N^ M), but the exponential map is no longer globally defined 
on this space since it's very easy for geodesies to intersect. See |EM| for details. 
We note that the metric dH) is invariant under the right-action of V^, (TV) and under 
the left-action of the isometry group of M. Hence even if M = TV, the metric is 
neither right- nor left-invariant on the open subset I?(M). 

li T] G C°° {N, M) is a diffeomorphism, then we can express any element U € 
TriC^{N, M) as [/ = uo 7], where w is a vector field on M. By the change of vari- 
ables formula for integrals, the Riemannian metric (j4]) becomes the more familiar 
expression 

{{uo-q,vori))r,= / p{q){u{q),v{q)) dp. 
Jm 

The geodesic equation = on C°°{N, M) becomes, in terms of u defined by 

(5) T]t{t,p) ^ u(t,7]{t,p)), 

the Burgers' equation (or pressureless Euler equation) ut + WuU = 0, while differ- 
entiating the density formula ([3]) leads to the continuity equation ([2]). 

To obtain the barotropic equation ([1]), we define (see Ebin |Eb| or Smolent- 
sev \Sm\ ) a potential energy function ^: C°°{N, AI) — > R which depends only on 
the density, of the form 

(6) »(,) = ^^„,Hrt*^X,"(34y)<i- 

Wc can then compute that the gradient of this function in the metric ([4]) is 
(grad$)^ = Q gradp(p)^ o jy, where p{p) = p^ ip' {p) , 

and hence Newton's equation — ^(gTad$),, can be written in the form ([T]) 

together with As pointed out by Smolentsev |Sm| . this system (like any 

conservative Newtonian system) can be rewritten as a geodesic equation with a 



^To prove theorems about geometry on infinite-dimensional manifolds rigorously, one prefers 
a Banach or Hilbert structure; however those theorems will not apply in this case regardless, so 
we lose nothing by staying in the C°° category. 
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modified metric — the Jacobi or Maupertuis metric. Generally speaking, if 7 sat- 
isfies Newton's equation = ^(grad on some manifold, then the energy 
E = ^|7(i)P + ^(7(^)) is constant in time, and all solutions with the same energy 
will be geodesies in the conformally equivalent Riemannian metric 

{u,v)^ = (£;-4>(7))(m,w)^. 

In this paper we will take a different approach which leads to a different formula 
for the Riemannian curvature. 

In recent years many authors have studied infinite-dimensional geodesic equa- 
tions which arise on groups of diffeomorphisms or related groups with right-invariant 
metrics. The geodesic equation in this case takes the form r](t) ~ dR,^(t-)u(t) and 
u{€) -I- ad*(j) u{£) = 0, and the second equation is the Euler-Arnold equation. See 
jAKj for examples. The closest relevant example in the present situation is the 
semidirect product V{S^^ k C°°(5^). The Euler-Arnold equation was computed in 
[S] to be 

ut + 3uu^. + ffx = 0, ft + fu^x + ufx = 0. 

Although this resembles the system ([T])-®, the extra factor of 3 makes it genuinely 
different, and this factor cannot be removed by a rescaling without changing the fiow 
equation rjt = u o rj. In addition, the higher-dimensional version as in [V| does not 
even superficially resemble the equation ([T]). More general Euler-Arnold equations 
on semidirect product groups were considered in j'HMR] , who also considered an 
alternative Lagrangian approach to the barotropic equations (where p is considered 
as an independent variable rather than derived from flow ?/) . Our approach uses a 
non-invariant metric which does not lead to an Euler-Arnold equation. 

3. The geometry of C°°{N,M x M) 

Now we define our (noninvariant) Riemannian metric. Our configuration space 
wiU be C°°(7V,M) x C°°(iV,M) = C°°{N,M x M). Tangent vectors to a point 
{t],F) are of the form (C/, 0) where U e T^C°°{N,M) and 4> e C°°(7V,M). If -q is 
a diffeomorphism, we can express (?7, (j)) = (u ° f ° v): where w is a vector field 
on M and / is a function on M . Let A : M+ IR+ be some smooth function, and 
define a metric on the product C°° {N, M x M) by the formula 

(7) {{{uovJ °v)Av°V,9°v))){n.F) = + u)] d//. 

Jm 

This is essentially a warped product of C°°{N, M) in the noninvariant metric (|3]) 
with the space of functions on N. 

3.1. Basic formulas. The following Lemma shows us how to differentiate func- 
tions on our configuration space. The proof is a computation which can be found 
in [Sm] or [Eb]. 

Lemma 3.1. Suppose W is a vector field on C°°{N,M x R) given by VF(r;.F) = 
(w o ri,ho T]) for some vector field w on M and function h: M — > M. Suppose 
$ : C°°{N, M xR) is a function of the form 



^{r],F)= / aip{p)dii— / a o Ja.c{r]) ip{l/ 3ac{r])) dv 
Jm Jn 
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for some functions a: M — > M and tp: — > K, where p is defined by ([3]). Then 
the derivative of <I> in the direction W is given at any point (77, F) by 

(8) W^ri,F)'^^ I P{'w,grad{a<p'{p)))dp.^ - div {pw)a<p' {p) dp. 

JM JM 

Now that we know how to differentiate functions that depend only on the density, 
we can use the Koszul formula to obtain the covariant derivative. For our purposes 
it is sufficient to compute in terms of vector fields of the form W(^ri,F) = {worj^horj) 
for some vector field w on M and function /i : Ai^ — !• R, as in Lemma [3. II 



Lemma 3.2. Let u and v be vector fields on M , and let f and g be real-valued 
functions on M. Define vector fields U and V on C°°{N,M x R) by U^ri,F) — 
{u or], f o 77) and V(,,_j;-) = {v o 77, g o 77). Then the covariant derivative VuV is given 
by (V;7y)(,,^i?) = Z or], where 

Z ^ {uig),Vuv) +Tp{iuJ),iv,g)) 

and 

Tp-. (x(M) X C°°(M,M))^ ^x(M) X C°°(Af,M) 
is the Christofjel map, a bilinear map satisfying 



(9) {{rp{{u, f),{v,g)),{w,h))) ^ I tp{p){hfdivv + hgdiYu-fgdiyw)dp, 
where p is given by Q and 



M 



(10) <p(p) = i(A(p)-pA'(p)). 

Proof. This is straightforward using the Koszul formula (see e.g., |dC| ). Lemma 
13. 1[ and the fact that the commutator of right-invariant fields is [J7, V^](^ i?) = 
077,(^(5) -«(/)) 077). □ 

Remark 3.3. Note that in general a weak Riemannian metric does not necessarily 
have any covariant derivative at all; what happens here is essentially the same sort 
of fortunate accident that occurs on the full diffeomorphism group T>{M) as in Ebin- 
Marsden |EM| . Here we are effectively working on C°°{N, M x M.), and using the 
fact that the covariant derivative constructed by jEM] on T>{M) works in the same 
way on C^{N,X) for any Riemannian manifold X. Of course the usual covariant 
derivative on C°°{N, A/xR) is somewhat different since that comes from a Cartesian 
product metric and our metric is a warped product; however the difference of two 
connections is an operator defined pointwisc (the Christoffcl symbol), and so we 
get our connection as long as we understand this difference. Lemma 13.21 basically 
just computes this difference. 

Although the formula (j9]) will be most convenient for our purposes, it is easy to 
see via integration by parts that we can write Tp more explicitly as 

(11) rp((uo?7, / 077), (7; 077,5 07;)) = (2077,^077) where 

j = Trlif + 9<i^vu) and z = - grad ((p{p) f g) . 
X{p) p 
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Notice that if A(p) = p, then f{p) ~ and hence Fp — as weU. In this case the 
geometry reduces to the geometry of M (as we will see later), and the covariant 
derivative is 

(12) {^uV){r,,F) = DR(^r,,F) {u{g),'^uv) = {u{g) o r], V„w o ij) 

when [/(,,, _F) ^ {uoT],f ori) and V(,,^f) = {v o -q, g o -q). 

This will be useful later. Having obtained the covariant derivative, we can start 
computing. Our first goal is to obtain the geodesic equation. 

Corollary 3.4. Suppose (^ri{t), F{t)) is a geodesic curve in C°°{N,M x M) in the 
Riemannian metric ([7]), with r]{t) a dijfeomorphism from N to M on some interval 
{—T,T). Define p[t) by ([3]). Let u[t) be the vector field satisfying ^ =uori, and 
let fit) be the function satisfying ^ = / ° ?/• Set q = \{p)f . Then u and q satisfy 
the equations 

(13) wt + V„u+- grad(4^) =0 

(14) + div (qu) = 0, 
where </? is given by pOp . 



Proof. If we define u and / by 77 (t) = u{t) o rj[t) and F{t) ~ f{t)° v(t) ^^'^ 
U{t) — (it(t), /(i)) , then the geodesic equation can be written as U{t) o rj{t) + 
(V(7(j)C/(t)) o rj{t) ~ 0. Using formula (|TT1) . this becomes 

^ + V„u + - grad ((^(p)./') = 

ft + uif) + ^^fdwu^O. 

Hp) 

Using it is then easy to see that q = Hp)f satisfies HH). □ 

Obviously (|14p is exactly the same differential equation as ([5]), which implies 
that if q = /9 at time zero, then q = p for all time. In this case, the equation for u 
satisfies 

.. + V„.+ igrad(4^).0. 

p ^ Kp) ' 

This is precisely equation ([T|) if we define 

pV(p) 



(15) p{p) = 



Kp? 



We have thus obtained the barotropic equations as a special case of the geodesic 
equation (fT3 |) -(fT4 |) . We summarize this as another corollary. 



Corollary 3.5. Suppose 70 = (r;o,-Po) is a point in C°°{N,M x R), and set po to 

be the density 0/770 given by Let Uq G T(^,-i^^ p^^-^C°°{N,AI x R) be the vector 
Uq = (uo o T^Q, j^jj ° Vq)- Then the geodesic ^{t) with initial position 70 and initial 
velocity Uq has tangent vector 'y'{t) = DRj(^t-)(^u{t), f(t)) , where u satisfies the 
barotropic evolution equation ([l} and f = for all time. 
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We thus obtain the equations of interest if we restrict the initial velocity in 
C°° (TV, M X R) to be something very specific in the function direction (but arbitrary 
in the diffeomorphism direction). This gives a distribution of allowable velocities 
in C°°{N,M X M). It is easy to see that this distribution is nonholonomic; in 
other words, there is no submanifold of C°°{N,M x R) for which all geodesies 
will correspond to barotropic flow. This limits the applicability of our standard 
geometric techniques, although we could if desired make sense of this situation in 
the context of the symmetric product] see for example |BLj and references therein 
for the general context. 

Remark 3.6. It is natural to ask whether the geodesic equations (fT3 |) -(fT4 |) have 
any meaning if our initial condition is not of the special type in Corollary (|3.5p . 
An interpretation is as follows. More general compressible fluid mechanics allows 
for the pressure p in ([1]) to depend on both the density p and the entropy s. In the 
absence of shocks, the entropy of any fluid particle is conserved, which implies the 
equation 



If instead of having q ^ p in ([M]) we have q ~ pC{s) for some function : M — > M, 
then the fact that pt + div (pu) = implies that s must satisfy (|16p . as expected. 
The corresponding pressure function would then be 



In other words, we can represent any compressible fluid as a certain family of 
geodesies in C°° {N, M x R) as long as the pressure is separable as a function of its 
arguments, and conversely every geodesic in C°°{N, M x R) is a compressible fluid 
flow for some choice of the separable pressure function. 

3.2. The sign of the curvature. In spite of Remark 13.61 we will still consider 
in the rest of this paper only barotropic flow. Since our primary motivation for 
considering the barotropic equations as geodesies is to understand stability in terms 
of curvature, we will compute the sectional curvature of this manifold. We are 
interested primarily in those sections where at least one of the vectors lies in our 
nonholonomic distribution, but we will first work out the general formula. 

Theorem 3.7. Let U and V he right-invariant vector fields on the manifold C°°{N, Mx 
R) given by f) = o r], f o rj) and f) = o r/,g o rf) for some vector fields 
u and V on M and functions f and g on M. Then the (unnormalized) sectional 
curvature of the metric ^ is given by 



(16) 



St + u{s) = 0. 



(17) 



PiP,s) 



Kp? 



(18) 




g div u] dp 



2fgQ(u,v)] dp 
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where the symmetric operator Q is defined by 

(19) Q{u, v) = div (Vijw) — u(divw) — (div w)(divti) 
and {R(u,v)v,u) is the curvature on M. 

Proof. This is a long but straightforward computation using Lemma 13.21 The fact 
that Q is symmetric fohows from the computation 

Q{u, v) — Q{v, u) = div [u, v] — u(divf ) + u(divu) ~ 0. 

□ 

The terms of are for the most part relatively simple, and thus we can 
determine the sign of the curvature quite easily, at least in special cases. 

Corollary 3.8. If N = M = and if 

(20) xip'(x) + ^4r^>0 forallx>0, 

X{x) 

then the sectional curvature (jlSp of C°°{N^M x M) is nonnegative in all sections. 

Proof. In one dimension the first and third terms of pSI) always vanish, so positivity 
of the second term is sufficient. □ 

Note that we cannot prove strict positivity: if /, g, u, and v are C°° functions 
with disjoint supports in 5^, we will have {{R{U, V)V, U)) = 0. 

Remark 3.9. For a polytropic fluid, where p{p) = Ap^ for some constants A and 
7 > 1, it is easy to check that the condition (|20p is equivalent to 7 < 3. This is 
true in many applications (e.g., for typical gases at room temperature, 7 = 1.4; 
later we will consider the special cases 7 = 3 in one dimension and 7 = 2 in two 
dimensions). 

In higher dimensions the formula gets more complicated, and ensuring positiv- 
ity of curvature is more difficult. We choose two examples for which it is easy to 
compute the Jacobi fields (and hence completely determine the linearized stability) 
in order to illustrate the relationship between positive curvature and stability. We 
emphasize that since the Riemannian exponential map is not smooth or even C^, 
we cannot rigorously prove any relationship between positive curvature and bound- 
edness of Jacobi fields using tools like the Ranch theorem, even for short time. 
Hence these computations should be viewed as intuitive guides rather than directly 
useful for stability analysis. 

Example 3.10. Suppose M = with fiat Riemannian metric and u is a velocity 
field of the form u = u {x) ^. Then for any pressure function p{p)^ the velocity 
field u is a solution of the steady compressible Euler equations 

(21) ViiU = — - gradp(/9), div (pii) = 

P 

with constant density p = 1. Writing U = {u,p/\{p)) and V = (y,p/\{p)), the 
curvature is given by 

(22) ((i?(t/, V)V, U)) = ^^(1) /^^ (div vf dx dy. 
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This works because VuU = for such a vefocity field, which forces p to be 
constant. Since Q(z, z) integrates to zero for any vefocity field z, the corresponding 
terms disappear. 

Example 3.11. Suppose M is the unit disc in M^. Let m = ^ be a rotational 
velocity field. Suppose that the pressure function is given by p{p) = for 
some constant c, so that A(p) = ^ and ^p{p) = 

Then u is a solution of the steady Euler equation (|2ip with density function 
pi''') = Eind the sectional curvature in directions U = (u, p/ X{p)) and V = 
(/cu, p/X{p)) for some constant k is given by 

iiRiU,V)V,U)) = -^^^^^- 
To see this, just notice that since f = g = p/X{p) in (fTBjl . it reduces to 

{{R{U,V)V,U)) = ^ [ p\divz)'dp+^ [ p^Qiz,z)dp, 

where z = v — u. Now div z = while Q{z, z) = — 2(fc — 1)^, and the rest is an easy 
computation. 



Remark 3.12. The geometric approach of Smolentsev jSmj . using the Jacobi met- 
ric, yields an alternative sectional curvature formula; let us compute its sign. For 
simplicity we will work in the one-dimensional case. In this case the configuration 
space is V^S^), and it is natural to compute the sectional curvature at an 77 € 2?(5'^) 
in a plane spanned hy U = u o rj and V — v o rj. We can suppose the fields are 
normalized so that {{U, V)) = Jg, p{x)u{x)v{x) dx ^ 0, and {{U, U)) = ((F, V)) = 1. 
After some simplifications, the formula of Smolentsev then yields that the sectional 
curvature K of T>(S^) in the Jacobi metric is given by 



K 



pp'{p){u" +v'^)dx 

2 



4(£;- $(7;))2 

3( /gi p'{p)p'v dxf + 3(/5i p'{p)p'u dxf - p'ipfp'^lpdx 



E - $(77) 

Clearly the fact that the first term involves (w'^ -I- u'^) implies that for typical u 
and V we obtain positivity. 

But we now ask whether this curvature can ever be negative. Now the parameter 
E is arbitrary, as long as £^ > <&(f?o): and thus the term 1/[E — ^{rf)] can be made 
arbitrarily large. In addition, for a given p and p, we can certainly choose u and v 
so that p'{p)p'v dx = Jgi p'{p)p'u dx = 0. Hence for any nonconstant density p, 
for values of E sufficiently close to $(77) there will be velocity fields u and v such 
that the sectional curvature is approximately 

We mention this only to emphasize that the Jacobi-mctric approach is funda- 
mentally different from our approach as far as curvature and stability predictions 
go, in spite of the fact that both generate the same geodesic equations. 
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4. Jacobi fields and linear stability of compressible motion 

In this section wc compute some examples of Jacobi fields explicitly, partly to 
illuminate the Lagrangian stability theory and partly to discuss the properties of 
the Riemannian exponential map. We work with the case p{p) = \p'^ on the circle 
(where the Euler equations ([I])-© reduce to the Burgers equation and can be 
solved fairly explicitly) and the case p{p) = ^^p^ in two dimensions. 

Note that although we have formulas for the covariant derivative ([TTjl and the 
curvature p^ . which in principle allow us to study the Jacobi equation directly, 
the formulas are complicated enough and there is enough cancellation in them to 
make them more trouble than they are worth. Instead we work directly with the 
linearizations of equations ([T])"©- 

Theorem 4.1. Consider a family C,{s, t) of geodesies in C°°{N, M x K), depending 
on a small parameter s, with tangent vectors satisfying the condition of Corollary 
\3.5[ with C(0, t) = {ri{t), F{t)) and ^(s, 0) = id for all s. Let u and p he the solutions 
of (dl)-®, related to r] through ^ and 

Let J{t) = f§(0,t) he the corresponding Jacohi field. Then J{t) = {j{t) o 
r](t), G{t)) for some vector field j (t) and function G{t), which satisfy cr = — div {pj) 
and G{t) — g{t) o r/(t) along with the linearized Lagrangian equations 

(23) g ^ 2if{p)a/X{pf + jp/ X{p) and — + [u, j] = v 
and the linearized Euler equations 

(24) ^ + div (cru) + div (pv) 

ot 

(25) — +V„i; + V„u + grad(/i'(p)CT) ==0, 
where h is the function defined hy h'{p) = p'{p)/ p. 

Proof. Write ({s, t) = (^(s, t),F{s, t)) , with u{s, t) and p(s, t) such that r]f ^uorj 
and 'poriJirf) = 1. Define a = ds\^_^ and v = ds\^_^. Differentiating ^ and 
^ with respect to s yields ([Ml) and ([25]). Differentiating dt'F = {p/X{p)) o rj with 
respect to s gives 

which reduces to the first part of (|23p . The second part of comes from differen- 
tiating the flow equation ([5]) and using the general formula = ff pointwise 
on AI (see e.g., do Carmo |dC| ). Finally the relationship = — div (pj) comes from 
differentiating equation □ 

4.1. The one-dimensional case. Consider the one-dimensional case M = 
(corresponding to periodic motion on R). We will assume p{p) = ^p'^; this is one 
of the cases that can be solved explicitly (see Courant-Friedrichs |CFj . Chap. Ill, 
Sect. 28), which allows us to write down all the Jacobi fields explicitly as well. 

Theorem 4.2. Suppose (rj^F) is a compressible geodesic as in Corollary \3.5\ with 
density p and velocity u on , with pressure function given hy p(p) = p^/3. Con- 
sider a Jacohi field along {rj, F) as in Theorem \4.1\ with initial condition J(0) ~ 
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and J'(0) ~ (wo,0). Then as long as the solution exists we have 

(26) IIjWIIl- <i||wo||L~. 

Proof. In one dimension with p{p) = /0'^/3, equations ([T]) and ^ become ut + uu^ + 
PPx ~ and pt + upx + pux = 0. We thus find that the functions 0+ = u + p and 
a_ = u — p both satisfy Burgers' equation: at + aa^; = 0. Let ^± be the flow of 
a±\ then of course we have 9(^^± = 0, the solution of which is obviously £,{t,x) = 
X + taQ{x). Hence a± satisfy the implicit equations ao{x) = a; + tao(a;)) . If x± 
denotes the spatial inverse of then we have a±{t,x) = ao(x±(^i2;)). 
The Jacobi equations in this case are 

(27) [pj)x = -o-, jt + ujx - jux = V, 

(28) at + {cru + pv)x = 0, vt + [uv + ap)x = 0. 

Define /?+ = v + a and /3_ = w — cr. Then the functions f3± satisfy the linearized 
Burgers equation 

A + + ^ 0. 
It is easy to see that the solution is 

I3±{t,x) = !'°},'^^]',^\ = vo{x{t,x))xx{t,x), 

since a{Q,x) = implies /3±(0,a;) = vo{x). Solving for cr and v, equations ^7} 
imply 



(29) pit,x)jit,x) = - / voiy)dy. 



X-it.x) 



X+{t,x) 



Now since the inverse flows x± satisfy x±(i, x) — x — ta±{t,x), formula ([^^ 
yields 

p{t,x)\j{t,x)\ < i sup\vo{y)\\x-{t,x) - x+{t,x)\ 
2 yesi 

< ^l|t'olU-|a+(i,a;) - a-{t,x)\ 
= tp(<,a;)||wo||L~, 

which implies (|26l). □ 



Recall that the curvature is nonnegative by Corollarv l3.8l but not strictly positive. 
In fact although our 7 = 3 is the critical case in ([20|) . which makes all terms but 
the last in (fT8|) vanish, the last term is generally positive, since we cannot expect 
the function / = p/A(p) to coincide with the function g satisfying Hence the 

curvature is nonnegative but sometimes positive, and the linear growth of Jacobi 
fields aligns with our intuition and is the best we can expect. This is a sort of weak 
instability: polynomial but not exponential growth of the Lagrangian perturbations. 

We can also use this explicit solution to establish the lack of smoothness of 
the Riemannian exponential map. Our technique involves analyzing the conjugate 
points along a particular geodesic, as in |M2| . jCKj . and |P4j . The idea is that al- 
though the exponential map is continuous in various function spaces (see Kato [K]). 
it cannot be in any Banach space: if it were, then the inverse function theorem 
and the fact that its derivative at is the identity map would imply the existence 
of a small interval on the geodesic in which no two points are conjugate (see e.g.. 
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do Carmo |dCj ). Hence we will prove the exponential map cannot be by finding 
a particular geodesic 7 such that j{tn) is conjugate to 7(0) for a sequence of times 
tn converging to 0. 

Theorem 4.3. Let H'{S^,S^ x R) denote the closure ofC°°{S\S^ x K) in the 
Sobolev topology, and consider the weak Riemannian metric with X{p) = 3/p. 
Then the Riemannian exponential map cannot be for any s > 3/2. 

Proof. By ([10]) and (HS]), this choice of A gives p{p) = as in Theorem l4?2l A 
particular solution of the equaitons is p(t,x) = u{t,x) = 1 for all t and x. Let 
7 = (ri,F) be the corresponding geodesic as in CoroUarv 13. 5( to prove that 7(0) 
and 7(r) are conjugate, we need to find a Jacobi field with J(0) = J{T) = 0. By 
theorem 14.11 we can write J{t) = (7 o rj, G) where j satisfies (P7)) and ^ ~ g o rj, 
and we can compute that g{t,x) = ^a{t,x). 

For this solution we have (as in the proof of Theorem 14. 2 1) that a+{t, x) = 2 and 
a-(t,x) = for all t and x, so that the inverse flows are x+i'^jx) = x — 2t and 
X- {t, x) = X. Thus we obtain 

a{t,x) = ^{yQ{x ~ 2t) — Vi^{x)) and v{t,x) — ^{yQ{x — 2t) + vq{x)Y 

Formula (HH) yields 

1 r 

3{t-,x) = 7; My)dy 

^ Jx-2t 

and since the flow of u is T]{t, x) ^ x + t, we obtain 

G{t,x)= 1 g{T,ri{T,x)) dr ^ ^ j {vo{x - t) ~ vo{x + t)) dr. 
Jo 'J Jo 

Now suppose vq{x) = cosn.T for some integer n. Then we have 

sin nt cos n( X — t) 

j{t,x) = ^ 

n 

G{t, x) = — sin (nx) sin^ {nt/2), 

and we see that there are conjugate points at T = for any positive integer m. 
Thus there is no neighborhood of on which the derivative of the exponential map 
is an invcrtible linear map. □ 

4.2. The two-dimensional case. We can do everything fairly explicitly in one 
dimension, at least assuming a simple pressure function. In two or more dimensions, 
explicit nonsteady solutions are much harder to come by. Hence we will just work 
out the Jacobi fields and Lagrangian stability for the two special steady solutions 
already studied in Section [321 via Proposition 13.101 and Example 13.111 

Theorem 4.4. Suppose M = T^. Suppose u = lo , for some constant lo, is a 
steady solution of the Euler equation (j2ip . with constant uniform density p = 1. Let 
c be the speed of sound, defined by c^ ~ for an arbitrary function p{p). Let 

7 = (?/, F) be the geodesic as in Corollary \3.5l Let J ^ [j o rj, G) be a Jacobi field 
along 7 as in Theorem \4.1\ with initial conditions j{0,x) = and jt{0,x) = vq{x) 
for some velocity field vq . Then j is bounded in time if and only if Vq is a gradient; 
otherwise j grows linearly in time at every point. 
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Proof. Equations (|23 l) -([25 |) reduce in this case to 

{dt + ujdy)j = V, [dt + ujdy)a + div v = 0, [dt + ujdy)v + grad a = 0. 
Ehminating a we get 

{dt + ujdy )^v ~ (? grad div 
with initial condition satisfying 

(30) i;t(0,a;,y) +a;i;y(0,a:,?/) = 0. 

It is easy to write down solutions: we expand v using the Hodge decomposition as 
V = grad / + where / is some mean-zero function and w is some divergence-free 
vector field which both depend on time. Then / and w satisfy the equations 

{dt + i^dyff = c^A/ and {dt + ujdyfw = 0. 

The solution satisfying (|30p is clearly 

f{t, X, ?/) = ^ an cos {cXnt)(j>n {x,y - Ldt) 
n 

w{t,x,y) = z{x,y-ujt), 

where A(/)„ = —\^^(j)n, the numbers a„ are some constants, and z is some divergence- 
free field. Integrating once to solve for j we obtain 

j{t,x,y) = > grad0„(a;,?/ - ujt) + tz{x,y - ujt). 

n 

Every component of the gradient part is bounded in time, while any nontrivial 
divergence-free part grows in time. □ 

Note that we obtain the behavior of Jacobi fields regardless of the curvature: in 
Proposition 13.1(11 we showed the curvature was nonnegative if and only if Lp'{\) + 
ip{\Y / \{1) > 0, which e.g., happens for pressure functions p{p) = Ap"^ if and only 
if 1 < 7 < 3. On the other hand even if this quantity is negative, it does not change 
the Jacobi fields in any way. 

Now we consider our final example: a rigidly rotating disc in the plane, as in 
Example 13.111 We assume the steady solution takes the form u ~ u) for some 
constant w, and that the pressure function is given by p{p) ~ so that the 

Euler equation ([T]) is ut -t- V„u = — c^grad/o for some constant c (the speed of 
sound). In this case we must have p'{r) = w^r/c^. Now the stability analysis 
depends on the boundary conditions we use. For our computations it is simpler to 
work with the free-boundary case (without surface tension), in which the density 
is specified on the boundary by p = po for some constant. If instead we worked 
with a fixed boundary (e.g., the fluid in a solid container), then the density would 
be unspecified at the boundary but the velocity would be required to be tangent 
to the boundary. 

We may assume by rescaling units that the boundary is at r = 1 and that po ^ 1. 
Then p{l) = 1 and p'{r) = up-rji? imply that p{r) = p^ - ^ + Since the 

2 

density must always be positive, we must have po > ^ for this to make sense. 
Hence the velocity may exceed the speed of sound but not more than by a factor of 
•\/2. Although in this supersonic region the evolution equations are not hyperbolic, 
and hence the standard existence theory breaks down, there is no apparent problem 
with the linearized equations in Theorem 14. II 
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Now we analyze the linear stability of the uniformly rotating fluid. Although 
stability can be analyzed for such fluids (and in much greater generality; see e.g., 
[HMRWj ) ■ we are interested in obtaining the explicit timc-depcndcncc of solutions 
of the linearized equations in order to study growth of the Jacobi fields as in The- 
orem 14.11 What is interesting about this computation is that we can express all 
the perturbations in terms of discrete Fourier modes: in general the fact that the 
corresponding operators are not self-adjoint means that wc should expect a con- 
tinuous spectrum as well as a discrete spectrum, but here the spectrum is purely 
discrete. The other reason this is interesting is that we don't require an ansatz 
for the growth of a perturbation; rather, the dependence is a consequence of the 
equations. 

The only boundary condition for the linearized equations is that cr = on the 
boundary, a result of p being prescribed. We note that Beyer and Giinther jBGj 
analyzed the Jacobi equations in this situation, in greater generality but less detail. 

Theorem 4.5. Suppose M is the unit disc in R^, with u = uj a rigid rotation 
for some constant ui. Further suppose that the barotropic fluid is described by the 
pressure function p(p) = (?p^ jl for some constant c. We assume the boundary 
conditions are such that the the density is a prescribed constant but the velocity is 
unconstrained. 

Then all solutions of the linearized equations (|24|) -([25]) are bounded in time, and 
we can write the solutions explicitly as sums of functions which look like Q(t, r, 9) ~ 
g^ytq(^j.'^gm{e-uit) j-^^ some real y, some n G Z, and some function q. 

Proof. The linearized equations are 

(31) at + div {pv + an) = 0, 

(32) wt + V„w + Vt,M = -c^gradcr, 

with boundary condition (7 = when r — \. To write these equations more ex- 
plicitly, we decompose the vector field pv into its gradient and divergence-free part 
(using the Hodge decomposition); more specifically 

(33) pv = grad / + sgradg = ^ ((/'^ + 9elr) dr + {fe/r^ - 9r/r) de^ . 

Without loss of generality we can assume that g vanishes on the boundary, although 
we cannot say anything yet about / on the boundary. 
We easily compute that 

VyU ~ — gradg sgrad/. 

r r 

Using the explicit formulas for p, u, and v, we can rewrite pip as 



(34) CT* + A/ + ^<Je = 0. 
Taking the divergence and curl, we obtain 

(35) CTt + were + F = 

(36) Ft + LoFg + 2uG + c^Acr = 

(37) Gt + i^Ge - 2ujF + uj'^ag = 0, 

where Act = div (p grader), F = Af, and G = Ag. By our free-boundary assump- 



tion, (T vanishes on the boundary, and compatibility requires that F and G have 
the same boundary conditions. 
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Since the hermitian operator A commutes with the antihermitian operator do, 
we can expand our functions in a mutual eigenbasis Cfcn('')e*"^, where A(<^fc„e™^) = 
— Afc„Cfcn('')6™^, for some functions (kn defined for k gN and n gN. By the usual 
Sturm-Liouville theory, we see that Afc„ — >■ oo as fc — >■ oo for any fixed n. 

Then we can write 

oo 

F(^,r,0) = ^^F,„(^)a„(r)e™^ 

k=l nSZ 

and similarly for G and a. Hence the coefficients satisfy the ordinary differential 
system 

dcTkn 



dt 



+ Fknit)=0, 



dFi 

(38) + 2ujGknit) - c'Afe„(Tfe„(t) - 0, 

- 2wFfc„(t) + mu^aknit) - 0. 

dt 

This is a constant-coefficient system, and its characteristic polynomial is + 
{c^Xkn + 4a;^)z + 2inuj^ = 0. and we want to show that it has three distinct 
imaginary roots in order to guarantee that a, F, and G are all bounded in time. 
Writing z = iy, we obtain 

(39) — Spy — 2q ^ 0, where p = {c^Xkn + 4a;'^)/3 and q = nuj^. 

Such a cubic has three distinct real roots if and only if p > and < p^. So 
we want to show that c^Xkn > uj'^{3n'^^^ — 4) for every n e Z and every k <E N. 
Obviously it is enough to show that c^Ai„ > a;^(3n^/'^ — 4). 

To obtain this bound on the smallest eigenvalue, we use the Rayleigh minimum 
principle. Let T be the space of smooth functions / : [0, 1] M such that /(I) = 0; 
then we can compute that 

. ^ J^\f'{rydr + n^J^'f{r)yrdr 

Xin > a mt — 



f^J" /o r/(r)2 dr 



. ^ J^r^f{r)^dr + n^J^rf{rrdr 
mi — i , 



where a = — ^ and b = The first term is exactly the Rayleigh quotient 
for the Bessel operator, and hence it is minimized when /(r) = Jn(c„r) with c„ 
the first positive root of the Bessel function J„; the minimum value is then c^. 
We claim the second quotient takes the minimum value n'^ + 1: to see this, write 
h{r) = rf{r), so that h{0) = h{l) = 0, and 

J^r^fir)'dr ^^^J^'rh'irrdr 



!lrf{rYdr ^KrY dr 



and the infimum of this last term is zero, as can be seen by computing it for 
functions h(r) ~ . ~ for a > 0; we obtain ^, so we can make it as close as 

desired to zero. As a result we have Ai„ > ac^ + h{v? + 1) > ^("-^ + !)■ The fact 
that c^Ain > ijj^['ir?l'^ — 4) now follows from the fact that + 7 > 671^/'^ for every 
integer n. 
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Hence all solutions of the system (p8| oscillate in time; specifically we can write 
cTknit) = crifc„e*^""* + cr2fcne*^2''"* + CT3fc„e'^="'"* for some distinct reals yikn, and 
similarly for Fkn and Gkn- Having obtained such a solution for F and G, we can 
then find / — A~^F and g = A^^G. Here g is assumed to vanish on the boundary, 
and we can obtain the boundary values for / using ((3T|) -(|32 ]) . □ 

It is now trivial to figure out how Jacobi fields grow, using (|23p . 

Corollary 4.6. Let 7 = {i],F) be a geodesic in the space C°°{D^,D^ x M) for 
which the corresponding tangent vector satisfies the condition of Corollarv \3.5[ with 
velocity u and density p given as in Theorem \4.5\ Then Jacobi fields are bounded 
for all time if the initial perturbation vq satisfies J^^ curl (^p{r)vo{r,9)) d9 = for 
every r. 

Proof. From (|23p we see that jt + ioje = v, and we expressed w as a sum of compo- 
nents of the form v{t,r,9) = e*y*e'"*^^~"*-'g(r) for some real number y. Hence the 
corresponding component of the Jacobi field is 

j{t,r,e) = e''^^^-'^'^q{r) [ e'^ ds. 

Jo 

This will be bounded for all time if and only if y 7^ 0. From equation ((39|) . we see 
that y = iff n = 0. It is easy to check that we have nonzero n = components if 
and only if the condition of the theorem is satisfied. □ 

This is another illustration of the fact that compressible fiows tend to be more 
Lagrangian stable than incompressible flows. There are many Jacobi fields which 
arc bounded in time in the compressible case, while in the incompressible case the 
curvature along a rigid rotational flow vanishes identically (it satisfies both the 
nonnegativity condition of Misiolek jMlj and the nonpositivity condition of the 
author |P3j ) and thus all Jacobi flelds grow linearly in time. 
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